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Known results on two-dimensional quantum electrodynamics (QED2) have been used to study the 
dependence of functional renormalization group equations on renormalization schemes and approx- 
imations applied for its bosonized version. It is demonstrated that the singularity of flow equations 
can be avoided in the optimized and power-law schemes for the bosonized model and the drawback 
of renormalization on bosonization is shown: it is indicated that renormalization of QED2 possibly 
requires interaction terms corresponding to higher frequency modes of its bosonized version. 
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I. INTRODUCTION 

In low dimensions, bosonization rules enable one to 
reformulate fermionic and gauge models in terms of el- 
ementary scalar fields. For example, the bose form of 
the two-dimensional massive Thirring model [l[ is the 
massless sine-Gordon (SG) scalar theory Q. Single fla- 
vor QED2 with massive fermions can also be bosonized 
and the corresponding scalar theory is the massive sine- 
Gordon (MSG) model [3|— 15| . Furthermore, the multifla- 
vor QED2, and the two-dimensional multicolor quantum 
chromodynamics (QCD2) can also be rewritten as multi- 
component SG theories [Hal- Thus, corresponding bose 
models are usually SG-type theories. The critical be- 
havior of original fermionic and gauge theories and their 
bosonized versions can be studied by various methods 
such as renormalization group (RG) approaches. If crit- 
ical behaviors of fermionic and gauge models are known, 
then bosonization transformations can be used to con- 
sider the dependence of methods suitable for the critical 
behavior of SG-type models on the approximations used. 

The goal of this paper is to study the dependence of 
functional RG equations obtained for the MSG model 
on renormalization schemes and applied approximations. 
The phase structure of the corresponding fermionic the- 
ory, the single flavor QED2 with massive fermions has al- 
ready been mapped out by density matrix RG approach 
Q which is used here to optimize the scheme-dependence 
for the MSG model. The drawback of the RG study of the 
MSG model on the bosonization transformations is also 
discussed, namely it is indicated that renormalization of 
QED2 possibly requires interaction terms corresponding 
to higher frequency modes of the MSG model. 

The paper is organized as follows. In Sec. [TIJ some 
aspects of bosonization are discussed. Functional RG 
equations are obtained for the MSG model in the sec- 
ond order of the gradient expansion in Sec. IIII1 Results 
of perturbative RG are summarized briefly in Sec. IIVI 
The nonperturbative RG study is given in Sec. [V] for the 
single-frequency, and in Sec. I VII for the multi-frequency 
MSG model in local potential approximation (LPA). In 
Sec. I VIII the RG flow is determined beyond LPA. Sec- 
tion |VIlT] serves as the summary. 



II. BOSONIZATION 

The mapping of quantum field theories of interacting 
fermions onto an equivalent theory of interacting bosons 
called bosonization is well-established in the context of 
1 + 1-dimensional theories. The well-known example is 
the massive Thirring model which is a theory of a 
single Dirac field if) determined by the Lagrangian density 

^Thirring = ^{vf^n ~ m)lp ~ -gfjfi (!) 

where j M = ipj^ip, m is the mass, and g is the coupling. 
This can be mapped onto the SG scalar field theory de- 
scribed by the Lagrangian density Q 
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where <p is a one-component scalar field and the iden- 
tifications Att/P 2 = 1 + g/ir, -P/{2n)e^d u (p = 3^ and 
u cos(/3tp) = —mipip, are made between the parameters 
of the two models. (Conventions and the definition for 
appropriate normal ordering are given in |2j.) The La- 
grangian of QED2 with a massive Dirac fermion which is 
also called the massive Schwinger model reads as Q 
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where F^ v — d^A^ — duA^. Using bosonization technique 
the fermionic theory (j3|) can be mapped onto an equiva- 
lent Bose form Q which is considered as the specific form 
of the MSG model 0, 0, Q whose Lagrangian density is 
written as 



£msg = ^(9^Y 
1 = 4tt, M 2 = 1 



+ - M > + u cos{(3ip) (4) 

with (3 2 = 4tt, M 2 = e 2 /7r, u = em exp (7)/(2tt 3 / 2 ) 
where 7 = 0.5774 is the Euler's constant and the vac- 
uum angle parameter has to be chosen as = ±7r for 
u > and Q = for u < 3- The MSG model has two 
phases. The Ising-type phase transition 0] is controlled 
by the dimensionless quantity u/M 2 related to the criti- 
cal ratio (m/e) c of QED2 which separates the confining 
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and the half-asymptotic phases of the fermionic model. 
The critical ratio (m/e) c = 0.31 — 0.33 has been calcu- 
lated by the density matrix RG method for the fermionic 
model which implies Q 



VhPJc = \TJ 



exp (7) 
2^ 



0.156- 0.168. 



(5) 



If one assumes a quartic self-interaction among the mas- 
sive Dirac fermions of QED2 by adding a Thirring type 
term to the Lagrangian ([3]) then one arrives at the mas- 
sive Schwinger-Thirring model which reads as 
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where F^ u = d^A v — d u A^ and — ipj^ip. It has 
been argued that by using bosonization technique the 
fermionic theory © can be mapped onto the MSG model 
if 4ir/f3 2 = l+g/(2w) and the Fourier amplitude is related 
to the fermion mass (u ~ m) and M 2 = e 2 /(it + g/2). 
Let us note that the bosonization of two-dimensional 
gauge and fermionic models (special attention on the 
Schwinger-Thirring model) has been the subject of in- 
tense study [9|. 



III. FUNCTIONAL RG METHOD FOR THE 
MSG MODEL 

In this section, we derive functional RG equations for 
the MSG model. Recently, the complete phase structure 
of the SG model ([2]) has been mapped out by extend- 
ing the functional RG analysis [Ty| beyond LPA fill ]. 
Here we use the same RG approach for the MSG model 
(fj|. Namely, the effective average action functional RG 
method [12Hl4j where the evolution equation reads as 



kd k T k = -Tr[(R k + r;')- 1 kd k Rk] 



(7) 



with the notation ' = d/dip, and the trace Tr stands for 
the integration over all momenta. As exact RG Eqs. ([7J 
are functional equations they are handled by truncations. 
Truncated RG flows depend on the choice of the regulator 
function Rk, i.e. on the renormalization scheme. In or- 
der to optimize the scheme-dependence various strategies 
have been worked out. For example, a general optimiza- 
tion procedure was proposed to increase the convergence 
of the truncated flow [Hj], and successfully applied in 
many cases, e.g. in quantum gravity [16j or in low-energy 
QCD [17]. Although, the optimized regulator i?£ pt [TJ] 
does not support the gradient expansion beyond second 
order [18l . ll9l ] but an optimization criterion based on func- 
tional variation is proposed to handle this problem [2(| 
and has been used for the study of low-energy behavior 
of the Yang-Mills theory. Since the RG study of the SG 
model was done by using the power-law type regulator 
Rf° w [l3j , it is a natural choice that for the RG analysis 



of the MSG model we use also the power-law regulator 
and the previously mentioned optimized one 
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, = a (P - p 2 )Q(k 2 - p 2 ) (8) 



where b > 1, usually a = 1 and O(x) is the Heaviside 
step-function. Let us note that various types of regula- 
tor functions can be used (e.g. the exponential one ]12\ ) 
but here we focus on those ([5]) which provide us the pos- 
sibility to perform momentum integrals of RG equations 
analytically. 

Another problem related to truncations is the singular- 
ity of RG flows. The appearance of spinodal instability 
(SI), i.e. singularity in the RG flow could be the conse- 
quence of a too drastic truncation. In this work we show 
that the singularity of the RG flow obtained for the MSG 
model can be avoided for the optimized and power-law 
regulators ([8]), if the RG equation obtained in LPA is 
integrated out directly i.e. using appropriate approxima- 
tions the effective potential remains convex [21(. Let us 
note that the MSG model has already been investigated 
by using the truncated Fourier expansion of the poten- 
tial which was found to be a too drastic simplifications of 
the functional subspace since SI appeared in the RG flow 
and the Maxwell construction resulted in a degenerate 
effective potential scheme- independently 0. 

Equation ([7]) has been solved over the functional sub- 
space spanned by the ansatz for the MSG model (H|) 
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where the local potential contains a single Fourier mode 
V k (<p) = ^M 2 (k) ip 2 + u(k) cosfa), (10) 
and the following notations are introduced 
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via the rescaling of the field ip ~ > ip/fl in Q and z(k) 
stands for the field-independent wave-function renormal- 
ization. Although RG transformations generate higher 
harmonics, we use the simple ansatz (|10[) first since in the 
case of the SG model it was found to be an appropriate 
approximation [llj. Then Eq. (UJ leads to the evolution 
equations [llj 



dkV k 
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V k kd k R k 



r V k "' 2 I V 2 kd k R, 



\dp 2 dp 2 



dV k 
dp 2 
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with V k = l/(zp 2 + R k + 14") and T = ^ir)- 1 J 27T dip as 
the projection onto the field-independent subspace. The 
scale k covers the momentum interval from the UV cutoff 
A to zero. Inserting the ansatz (fit))) into Eqs. (|12p and 
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(fT3| the flow equations for the coupling constants are 
(similar RG equations obtained for the SG model in (Tpj ) 
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with P = zp 2 + M 2 + Pfc where 9feM = 0. In general, the 
momentum integrals have to be performed numerically, 
however, in some cases analytical results are available. 
Indeed, by using the power-law type regulator function 
with 6=1 (i.e. the Callan-Symanzik scheme), the mo- 
mentum integrals can be performed and the RG equa- 
tions reads as, 



(2 + kd k )u 
kdkZ 
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24tt [(1 + M 2 ) 2 -u 2 ]' 



(2 + kd k )M 2 = 0, (16) 
with dimensionless couplings u — k~ 2 u, M 2 = fc~ 2 M 2 . 



equivalences, the approximated RG equations (3.17) of 
23] can be written as 
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with «2 constant. Perturbative results (Eqs. (|18|) 
and ([17]) ) demonstrate that the Kosterlitz-Thouless- 
Berezinski [25j (KTB) type phase transition known to 
take place for the SG theory [HI, disappears in case of the 
MSG model. Indeed, one expects an Ising-type second 
order phase transition for the MSG scalar model which is 
controlled by the dimensionless ratio u/M 2 — u/M 2 . For 
example, according to Eq. (fTfl) , in the infrared (IR) limit 
(k 0) the mass term becomes large M 2 ~ k~ 2 which 
freezes out the evolution of z and the Fourier amplitude 
has a trivial tree-level scaling u ~ fc~ 2 ; consequently, the 
ratio u/M 2 tends to a constant value. 

However, the perturbative RG flow produces arbitrary 
large ratio depending on the initial conditions (there is no 
upper bound), hence no critical value can be determined. 
Consequently, the perturbative RG flow is not suitable 
for the prediction of the Ising-type phase transition of 
the MSG model. 



IV. PERTURBATIVE RG 

Let us first consider the massless limit (M — > 0) when 
the RG Eqs. (fTl))) reduce to those derived for the SG 
model [11]. The spontaneously broken phase of the SG 
model is known to be equivalent to the neutral sector of 
the massive Thirring model (JT]). Indeed, the renormal- 
ization of the massive Thirring model has already been 
discussed and it was demonstrated that the scaling of the 
fermion mass (see, e.g. !22j), m — mo(/i/ A)(~ 9 ^ 9+7T ^ is 
identical to the solution of the linearized form of the RG 
Eqs. (jTHJ) if kd k z = and M = arc assumed and the 
equivalences g/(g + 7r) = 1 — 1/(zAtt), k = fi, u = m/i/ir 
and ua = m A/7r are used. 

For nonvanishing mass (M 7^ 0), it is illustrative to 
compare the perturbative RG equations given in (23l . |24T ] 
to that obtained by the linearization of Eq. (fT6|) which 
reads as 



(2 + kd k )u = 
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(2 + kd k )M 2 = 0, 



(17) 



Let us relate the parameters of (23j /it, k and m 2 to that 
of defined in Eqs. © and (fTU)) . namely u = fi, z = 
1/(47t 2 k), M 2 = m 2 /(4-7TK) and k — A. By using these 



V. FUNCTIONAL RG STUDY OF THE 
SINGLE-FREQUENCY MSG MODEL 

Exact RG equations (with a single Fourier mode) show 
a different picture since one cannot obtain arbitrary large 
IR value for the ratio u/M 2 . Let us first consider the 
RG flow obtained by the optimized regulator function 



© with a = 1 in the LPA (i.e. z = 1//3 2 
which reads as 
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In the broken symmetric phase, the RG trajectories 
merge into a single trajectory in the deep IR region which 
is characterized by the critical ratio [u/M 2 ] c = 0.0625 
(for fi 2 = Ait) and serves as an upper bound (see Fig.Q]). 
The critical value obtained by Eq. (p~9|) is less than the ex- 
act result ([5]), therefore it requires further improvement. 
Let us try to improve it by the optimized regulator ([8jl 
with a ^ 1 which has the following form in LPA 
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(2 + kd k )M 2 = 0. 
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(20) 
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FIG. 1: Phase diagram of the MSG model for /3 2 = 4tt. 
RG trajectories are obtained by the integration of Eq. (|19|l . 
Since SI does not occur in the RG flow, the critical ratio of 
the MSG model can be determined, [u/M 2 ] c = 0.0625. The 
arrows indicate the direction of the flow. 



However, for a ^ 1 spinodal instability (SI) appears in 
the RG flow in the broken symmetric phase, i.e. RG 
equations become singular in the IR limit and the RG 
flow stops at some finite scale (see the dashed lines in 
Fig. [2}. Although RG trajectories start to converge into 
a single one in the broken phase the critical value of the 
single-frequency MSG model cannot be determined un- 
ambiguously. In other words, the convergence properties 
of the optimized RG is weakened for a^l. Let us try to 
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FIG. 2: Phase diagram of the MSG model for /T = 4tt. RG 
trajectories are obtained by the integration of Eq. (|20[) with 
a — 100. The dashed lines correspond to RG trajectories 
where SI occurs in the RG flow, thus the critical ratio of the 
MSG model cannot be obtained. 



use other types of RG equations, for example the power- 
law type RG with 6 = 1 in LPA which reads as 



(2 + kd k )u 



(21) 

The RG Eq. (|2"Tj) can be obtained by rescaling u — > u/z 
and M 2 -> M 2 /z in (JTSJ) and using the identifications 
(|lip with the assumption dkZ = 0. It is known that in 
the sharp limit, the optimized RG becomes identical to 
the power-law type RG with 6=1, i.e. Eq. (|20[) reduces 
to Eq. ([2~Tj) for a —> oo. Thus, SI is expected in case of 
Eq. (f2"Tj) . Indeed, the numerical solution of (|2"Tj) indicates 
the appearance of SI in the broken symmetric phase (see 
the dashed lines in Fig. [3]). 
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FIG. 3: Phase diagram of the MSG model for /3 2 = 4tt. RG 
trajectories are obtained by the integration of Eq. (|21|l . The 
dashed lines correspond to RG trajectories where SI occurs 
in the RG flow, thus the critical ratio of the single- frequency 
MSG model cannot be obtained. 

It is also illustrative to compare the IR values of the 
ratio u/M 2 at the scale of SI given by the integration 
of optimized RG with various values for the parameter 
a using the same UV initial condition (see Fig.|4]). This 
demonstrates that the best estimate for the critical ratio 
of the single- frequency MSG model, in the framework of 
the optimized RG can be achieved for a = 1. Our find- 
ings are consistent to the feature of the optimized RG 
namely that it increases the convergence properties of 
the truncated flow. For example, similar result is shown 
in Fig. 12 of (2|| in the framework of the O(N) symmet- 
ric scalar theory in d = 3 dimensions which has been the 
subject of intense study on scheme-dependence (see e.g. 
27]). Let us note that the optimized RG with a = 1 
produces reliable results for the single-frequency MSG 
model in LPA as opposed to the Wilson-Polchinski RG 
;28] which was found to be inappropriate Q for the de- 
termination of the phase structure of the MSG model. 
This is a counterexample for the statement, namely, that 
the Wilson-Polchinski RG and the optimized one always 
provide us with the same critical behavior in LPA [18| . 
The mapping between the two latter RG methods works 
only if the potential is nondegenerate but this is not true 
for the broken symmetric phase of the MSG model. 

Since Eq. (|19|) has no singular behavior, the appear- 
ance of SI is expected to be the consequence of an inap- 
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a=l optimized RG in LPA 
a^l optimized RG in LPA 
b=l power-law RG in LPA 




0.4 0.6 

x=a/(a+l) 

FIG. 4: This figure shows how the IR value of the ratio u/M 2 
obtained by the integration of RG Eqs. (|19[) and (|20|) depends 
on the parameter a of the regulator function. The same ini- 
tial condition has been used for the numerical integration, 
(m(A) = 10" 5 ,M 2 (A) = 1(T 9 ). SI has occured for a 1, 
thus the RG flow stops at some finite scale where the ratio 
has been read off and plotted. It is possible to avoid SI but 
only for a = 1. For a — > oo (i.e. x — >• 1) the RG Eq. ((20)) be- 
comes identical to that was obtained by the power-law type 
regulator with 6 = 1 (|21[) ; consequently in this case the IR 
values of the ratio coincide. 



propriate approximation, e.g. too drastic simplification 
of the functional subspace. Let us also note that SG- 
type models undergo an infinite order (or topological) 
phase transition and it was shown fH\ that the study of 
a single-frequency model is sufficient to recover the KTB- 
type critical properties (higher harmonics were found to 
be irrelevant). However, the MSG model has an Ising- 
type second order phase transition, hence there is no rea- 
son to focus on the study of a single-frequency model. 
Consequently, in order to obtain reliable results and to 
avoid SI in case of the MSG model one has to incorporate 
higher harmonics generated by RG Eqs. (fl"12j). (fTBf which 
is discussed in the next section. 



VI. FUNCTIONAL RG STUDY OF THE 
MULTI-FREQUENCY MSG MODEL 

We now turn to the discussion of the MSG model 
including higher harmonics. Let us consider the RG 
Eq. (fl"2")l in LPA (i.e. d k z = 0) using the ansatz 



1 °° 
V k (<p) = -M 2 y 2 + ^2u n (k)cos(n(3ip). 



(22) 



71=1 



There are two ways to determine the IR scaling of the 
MSG model (i) either the partial differential Eq. (O 
has to be solved directly by e.g. a computer algebraic 
program using the initial condition Q (higher harmon- 
ics are generated by RG equations), (ii) or one can find 
the solution of ordinary differential equations given for 



the coupling constants which are obtained by inserting 
the ansatz (f2"2"|) into Eq. (IT2l . Let us first discuss the lat- 
ter case when it is unavoidable to implement a further 
approximation besides the LPA, namely the truncation 
of the Fourier expansion of the potential. In this case, 
the RG flow on the trajectories started at /3 2 = 4tt in 
the broken symmetric phase develops SI at some finite 
scale fcsi > M [8]. The way to go beyond fcgi is the 
Maxwell construction (i.e. the tree-level blocking rela- 
tion |29() which represents too strong constraint on the 
RG flow and results in a scheme-independent infrared 
value for the critical ratio [u/M 2 ] c = 0.159 M. How- 
ever, SI occurs in the RG flow as an artifact [aUll due 
to the truncated Fourier-expansion applied to the almost 
degenerate blocked action of the MSG model, thus, one 
has to solve directly the RG Eq. (|12[) in order to obtain 
reliable results. 

Indeed, the truncation of the expansion of the blocked 
potential in a series of base functions may become unreli- 
able when the blocked action becomes almost degenerate, 
i.e. when k 2 + V k n approaches zero [3(| ■ This motivates 
a direct numerical solution of the RG Eq. (|12l) for the 
blocked potential which avoids any assumption on the 
functional subspace where the solution is sought for and 
any truncated series expansion in some base functions. In 
this case, the appearance of SI is avoided and the critical 
ratio can be determined scheme-dependently. For exam- 
ple, Eq. (fl"2|) for the power-law type regulator with b = 1 
reads as 
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of which direct integration results in [u/M 2 ] c = 0.148 
(see the solid lines in Fig. [5J. Moreover, one obtains a 



0.25 - 



!£S 0.15 




0.05 



FIG. 5: Phase diagram of the MSG model for /3 2 = 4tt. RG 
trajectories are obtained by the direct numerical solution of 
Eq. (I23p (solid lines) incorporating the effect of higher har- 



better result [u/M 2 ] c — 0.159 for the optimized regu- 
lator and also for the power-law type one with 6 = 2. 
Consequently, in case of the direct integration of RG 
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equations (when the effect of higher harmonics are in- 
corporated), (i) SI can be avoided in the RG flow, (ii) 
the critical ratio is found to be scheme-dependent, (iii) 
the RG schemes can be optimized using the known re- 
sult (|5|) obtained for QED2. Let us note that the best 
result for the critical ratio (the maximum which can be 
achieved by any RG scheme in LPA) is obtained by reg- 
ulators which have good convergence properties, e.g. in 
case of the optimized RG one has to choose a = 1 and 
for the power-law RG the good choice is 6 > 1. Those 
regulators which provide for RG equations poor conver- 
gence properties like the optimized regulator with a =/= 1 
and the power-law type regulator with 6=1 are not suit- 
able to recover the best result for the critical ratio. For 
comparison, see Table HI 



opt. single-fr. 


opt. multi-fr. 


b=l multi-fr. 


b=2 multi-fr. 


0.062 


0.159 


0.148 


0.159 



TABLE I: Critical ratio obtained by various RG meth- 
ods. 




FIG. 6: RG trajectories are obtained by the integration of 
Eq. (|16[) . Dashed lines correspond to RG trajectories where 
SI occurs, thus the critical ratio r£ =1 cannot be determined 
unambiguously. Vertical line shows the dimensionful mass 
scale which remains unchanged under RG transformations. 



where c(6) 
the ratio 



6/(6 - l) 1 - 1 ^. Therefore, in the IR limit 



(25) 



VII. FUNCTIONAL RG BEYOND THE LPA 

Can the appearance of SI be avoided by the inclusion 
of the wave -unction renormalization? What is the IR 
value for the frequency (i.e. f3 2 = l/z) in the broken 
symmetric phase of the MSG model if the wave-function 
renormalization is included? How does the wave-function 
renormalization affect bosonization? In order to clarify 
these issues, let us go beyond the LPA and solve the flow 
Eqs. ([T4]) and ([f"5j) obtained for the single-frequency MSG 
model where the wave-function renormalization z is kept 
scale-dependent. Since the optimized regulator does not 
support the derivative expansion beyond second order in 
this section we focus on the power-law RG. 

Let us first consider the Callan-Symanzik scheme (i.e. 
the power-law RG with b = 1 ) where the flow Eqs. (fT4|) 
and (fT"S"|) reduce to Eq. (fT6|) . As is seen in the inset of 
Fig. [51 SI appears in the RG flow in the broken symmetric 
phase, i.e. the RG equation becomes singular in the IR 
limit and the flow stops at some finite momentum scale 
(see the dashed lines in the inset of Fig. [5]) . 

Since the convergence properties of the power-law RG 
are increased for b > 1, Eqs. ([14]) and (fl"5|) are solved 
numerically for b = 2 (see Fig. [7]). Independently of the 
actual value of b, the potential was found to become de- 
generate in the broken symmetric phase and the RG flow 
is determined by the degeneracy condition (similar re- 
sults were obtained for the SG model in [111') 



c(6)^ 



-i/b 



u + M = 
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tends to one (i.e. rf,(fc — > 0) = f) in the broken phase 
(see the dashed lines in Fig. [7] for b — 2). Therefore, the 
ratio becomes universal in the broken phase. In the sym- 
metric phase it tends to a constant IR value depending 
on the initial conditions. The critical ratio rg(fc) which 
separates the phases of the single-frequency MSG model 
is represented by the thick solid line in Fig.[7J Similar 
results can be obtained for b = 1 (see the dashed lines in 
Fig. [6]), but due to the poor convergence properties of the 
Callan-Symanzik scheme, SI appears in the RG flow in 
the broken symmetric phase and the ratio cannot reach 
its universal value. Let us note that the RG flow always 
stops at a finite momentum scale in the broken phase 
independently of b but a better convergence is obtained 
for b > f . This indicates that a more accurate calcula- 
tion requires the inclusion of higher harmonics as it was 
demonstrated in LPA. 

In general, the single- frequency approximation is "im- 
proved" by the inclusion of the wave-function renormal- 
ization. For example, the critical exponent v of the MSG 
model can be obtained in the framework of the power- 
law RG with b > f if z(k) is kept scale-dependent. It 
is known [7| that the MSG model belongs to the two- 
dimensional Ising universality class, thus the correlation 
length is a power-law function of the reduced tempera- 
ture £ ~ t~ v with v — f . Indeed, if one defines the corre- 
lation length in the symmetric (disordered) phase by the 
constant IR values of the ratio, £ ~ [1 — rb(k — > 0)] -1 
and the reduced temperature is given by the initial UV 
(k = A) values, t = MA)" 1 - rg(A)-" 1 ]/rg(A)- 1 then 
v = 1 is obtained (see solid lines in Fig. [7] for 6 = 2). Let 
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b=2 power-law RG with z 




k/A 



10" 5 10" 4 10" 3 10" 2 10" 1 1 

k/A 

FIG. 7: RG trajectories are obtained by the numerical in- 
tegration of Eqs. (|14|l and (|15|) for the power-law regula- 
tor with 6 = 2. Dashed lines correspond to RG trajecto- 
ries in the broken symmetric phase of the single-frequency 
MSG model. Vertical line shows the dimensionful mass scale 
which remains unchanged under RG transformations. The 
inset shows the scaling of the wave-function renormalization 
in the two phases. 



us note that in following Ref. the correlation length 
can be defined as £ ~ (M — /c c ) _1 in the broken phase 
where k c represents the momentum scale at which the 
ratio (|25|) becomes constant during the RG flow. If the 
reduced temperature is t = [r&(A) — r{j(A)]/rg(A) then 
one obtains again the power-law behavior with v = 1. 
Consequently, the RG equations derived for the single- 
frequency MSG model beyond LPA are sufficient to indi- 
cate that the model undergoes a second order Ising-type 
phase transition (while it is known that the SG model 
has an infinite order, KTB-type phase transition [HI). 

Let us consider the RG evolution of the wave-function 
renormalization which is equivalent to the inverse fre- 
quency, i.e. z(k) = l//3 2 (fc). In the symmetric phase 
z(k) becomes a constant in the IR limit depending on 
the initial conditions (see the solid lines in the inset of 
Fig. [7]). In the broken symmetric phase, however z(k) 
runs into infinity for k — > (see dashed lines in the inset 
of Fig.[7j), i.e. it has a universal behavior in the broken 
phase thus /3(fc) tends to zero. Therefore, if one assumes 
that bosonization identifications between the parameters 
of the fermionic and the corresponding bosonic theory 
hold also for the blocked action then our result has a 
drawback on bosonization, namely it indicates the neces- 
sity to construct the fermionic counterpart of the MSG 
model for 1 ^ An. 

Finally, let us mention an open question related to the 
renormalization of the MSG model. By the inclusion 
of the wave-function renormalization it was possible to 
determine the critical exponent of the correlation func- 
tion which demonstrates that the MSG model belongs to 
the universality class of the two-dimensional Ising model 
i.e. the O(N) symmetric scalar theory for N = 1 and 



d = 2. It is known that higher order polynomial terms 
are needed in case of the O(l) model in order to obtain 
the critical behavior at the Wilson-Fisher fixed point in 
a reliable manner. Thus, it might indicate that the crit- 
ical behavior of the MSG model found at the IR limit 
is the consequence of a Wilson-Fisher type fixed point 
which could possibly appear in the MSG model if higher 
polynomials of the field are incorporated in the RG flow. 
Therefore, it is an interesting open question to consider 
the phase structure of the MSG model with the inclu- 
sion of higher order monomials of the field, although 
these terms are not generated by RG equations using 
the ansatz ((4]). 



VIII. CONCLUSIONS 

Known results on QED2 have been used to optimize 
RG schemes for its bosonized version, the MSG model 
and to consider how the results obtained by RG equa- 
tions depend on various approximations used. By the 
inclusion of the wave-function renormalization and the 
direct integration of RG equations derived for the MSG 
model, we went beyond the previously used approxima- 
tions. 

It was shown that the inclusion of higher harmon- 
ics and the direct integration of RG equations are both 
needed to avoid the appearance of singularity in the RG 
flow of the MSG model (for the optimized and power- 
law regulators) and to recover the critical ratio of QED2 . 
It is also demonstrated that the optimized RG predicts 
a reliable result for the single-frequency MSG model in 
LPA and it is known that the Wilson-Polchinski flow is 
not suitable to map out its phase structure @. Thus it 
shows that the two latter RG methods do not produce 
the same critical behavior in LPA if the blocked potential 
becomes degenerate which is the case for the MSG model 
in its broken phase. 

Moreover, as a result of optimization, the best result 
for the critical ratio of QED2 (the maximum which can 
be reached in LPA) is obtained by those regulators which 
have good convergence properties, such as the optimized 
RG with a = 1 and the power- law RG with b > 1. Regu- 
lators with poor convergence properties (e.g. optimized 
RG with a > 1 and power-law RG with 6=1) are found 
to be unable to recover the best result. 

If one assumes that the identifications between the pa- 
rameters of the fermionic and the corresponding bosonic 
theory holds also for the blocked action then results 
on the MSG model indicate that the renormalization of 
QED2 possibly requires interaction terms which corre- 
spond to higher frequency modes of its bosonized version. 
The renormalization of the wave-function in the RG flow 
of the MSG model has also a drawback on bosonization. 
Since in this case z = 1//3 2 is scale-dependent, it is a ne- 
cessity to construct the fermionic counterpart of the MSG 
model for f3 2 ^ Att. Indeed, if one assumes a quartic self- 
interaction among the massive Dirac fermions of QED2 
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by adding a Thirring-type term to the Lagrangian Q 
then one arrives at the massive Schwinger-Thirring model 
Q which was proposed as the corresponding fermionic 
theory of the MSG model for /3 2 f 4tt. 

Finally, let us note that the scenario discussed in this 
work can possibly be extended for many-body condensed 
matter systems [3l|. Another open question related to 
the present work is the direct comparison between flows 
of fermionic and bosonic models. For example, one can 
try to compare the functional RG study of the SG scalar 
theory and the fermionic Thirring model which can pos- 
sibly be achieved by the extension of the RG analysis of 
the three-dimensional Thirring model [32J to the two- 



dimensional one. Furthermore, bosonized versions of 
multiflavor QED2 and mutlicolor QCD2 are also SG-type 
models, hence they represent further examples where 
the drawback of RG results can be studied on bosoniza- 
tion (including wave-function renormalization as the fre- 
quency becomes scale-dependent). 
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